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ABSTRACT 

The energy spectra of TeV gamma-rays from blazars, after being corrected for inter- 
galatic absorption in the Extragalactic Background Light (EBL), appear unusually 
hard, a fact that poses challenges to the conventional models of particle acceleration 
in TeV blazars and/or to the EBL models. In this paper we show that the internal ab- 
sorption of gamma-rays caused by interactions with dense narrow-band radiation fields 
in the vicinity of compact gamma-ray production regions can lead to the formation of 
gamma-ray spectra of an almost arbitrary hardness. This allows significant relaxation 
of the current tight constraints on particle acceleration and radiation models, although 
at the expense of enhanced requirements to the available nonthermal energy budget. 
The latter, however, is not a critical issue, as long as it can be largely compensated by 
the Doppler boosting, assuming very large 30) Doppler factors of the relativistically 
moving gamma-ray production regions. The suggested scenario of formation of hard 
gamma-ray spectra predicts detectable synchrotron radiation of secondary electron- 
positron pairs which might require a revision of the current "standard paradigm" of 
spectral energy distributions of gamma-ray blazars. If the primary gamma-rays are of 
hadronic origin related to pp or p7 interactions, the "internal gamma-ray absorption" 
model predicts neutrino fluxes close to the detection threshold of the next generation 
high energy neutrino detectors. 
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1 INTRODUCTION 

The recent reports on detections of very high energy (VHE) 
gamma-rays from blazars with redshifts z 0.1 (for a re- 
view see e.g. Hinton (2007)) initiated renewed debates on 
the interpretation of TeV gamma-ray spectra of blazars, 
in particular in the context of the level of the diffuse ex- 
tragalactic background radiation at optical and infrared 
wavelengths, often called also as Extragalactic Background 
Light (EBL). Initially, the tight link between these two 
topics - TeV blazars and EBL - became a subject of hot 
discussions prompted by multi (up to 20) TeV gamma- 
rays detected from a nearby BL Lac object, Mkn 501 
l| Aharonian et al. 1999f ). and by the reports claiming de- 
tection of high fluxes of EBL at far infrared wavelengths 
(lHauser et al. 19981: ISdilegel et al. 19981 : lLagage et al. 19991 : 
iFinkbeiner et al. 2000f ). However, it was quickly recognised 
that these two claims hardly could be compatible within any 
standard model of TeV blazars (see, for a review, Aharonian 
(2001)). 

A distinct feature of extragalactic gamma-ray as- 



tronomy is that VHE gamma-rays emitted by dis- 
tant (^ 100 Mpc) objects arrive after significant ab- 
sorption caused by their i nteractions wi t h EBL via 
the process 77 — > e+e" (|Nikishov 19621 : Ijellev 19661 : 
iGould and Schreder 19671 ). The reconstructed, i.e. the 
absorption-corrected gamma-ray spectrum from a source at 
a redshift z, Jq{E) = J h s {E)e T( - E,z ^ depends on the flux and 
energy spectrum of EBL through the optical depth t(E,z). 
Thus, at energies where t(E,z) 1, the primary gamma- 
rays suffer strong spectral deformation. 

The EBL consists of two emission components produced 
by stars and partly absorbed/re-emitted by dust through- 
out the entire history of galaxy evolution. As a result, two 
distinct bumps are expected in the spectral energy distribu- 
tion (SED) of EBL at near infrared (NIR) and far infrared 
(FIR) wavelengths, with a mid-infrared (MIR) "valley" be- 
tween these two bumps (see e.g. Hauser and Dwek (2001)). 
Generally, for almost all EBL models, t(E) is a strong func- 
tion of energy below 1 LeV and above 10 TeV; between 1 
and 10 TeV the e nergy-dependence of t(E) is much weaker 
l| Aharonian 200 lb . Consequently, one should expect signif- 
icant distortion of the VHE spectra of blazars at energies 
below 1 TeV and above 10 TeV, provided that at these 
energies r > 1. One can re- formulate this statement in a 
different way. Namely, for a standard ("decent") intrinsic 
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gamma-ray spectrum, the observer should detect very soft 
(steep) spectra at energies below 1 TeV and above 10 TeV 
from objects for which r ^ 1 at corresponding energies. This 
condition is safely satisfied, given the constraints on the min- 
imum EBL flux imposed by galaxy counts, for blazars with 
redshifts z > 0.15 like 1ES 1101-232 and for nearby ob- 
jects with 2 ~ 0.03 like Mkn 501. Even though the detected 
gamma-ray spectra from both objects in the correspond- 
ing energy intervals are indeed quite steep with a photo n 
index ~ 3 (I Aharonian et al"""l 999; Ah aronian et al. 2006aT ). 
they appear not sufficiently steep to compensate the func- 
tion f(E) — e T ^ E \ and thus to prevent a robust conclusion 
that the intrinsic VHE gamma-ray spectra of these blazars 
are unusually hard. 

In the case of Mkn 501, the intrinsic spectrum has 
a "non-standard" shape with a possible pile-up above 10 
TeV which has bee n interpreted as a "IR backg round - TeV 
gamma-ray crisis" ijProtheroe and Mever 200d ) or a need to 
invoke dramatic assumptions like a violation of the Lorentz 
invariance (see e.g. Kifune (1999)). However, a more prag- 
matic view which presently dominates in both infrared and 
gamma-ray astronomical communities, treats this "crisis" 
as somewhat exaggerated, especially given the ambiguity 
of extraction of the truly diffuse extragalactic FIR compo- 
nent from the much higher backgrounds of local origin (see 
e.g. Hauser and Dwek (2001)). Nevertheless, the recently re- 
ported low limits on the EBL at mid infrared wavelengths 
from the Spitzer deep cosmological surveys appeared quite 
high, for example at 70 /xm the EBL flux should exceed 
> 7.1 ± 1.0 nW/m 2 s ()Dole et al. 20061 ). This implies that 
the problem is not yet over, and one may still face a chal- 
lenge with the interpretation of the energy spectra of Mkn 
501 and Mkn 421 in the multi-TeV energy domain. 

On the other hand, the recent detections of TeV 
gamma-rays from blazars with redshifts z ^ 0.15 renewed 
the potential problems and challenges for standard models of 
TeV blazars. This time the issue has a more solid experimen- 
tal background, because the gamma-ray spectra corrected 
for the intergalactic absorption appear very hard ("harder 
than should be" ) even for the minimum possible EBL fluxes 
at optical and NIR wavelengths. Namely, the HESS collab- 
oration reported, based on the detection of TeV gamma- 
rays from the BL Lac object 1ES 1101-232, that any sig- 
nificant deviation from the lower limits of EBL determined 
by the int egrated light of galaxies re solved by the Hubble 
telescope (IMadau and Pozzetti 20001 ) . would lead to very 
hard intrinsic gamma-ray spectru m with a slope character - 
ized by a photon index To ^ 1.5 ([Aharonian et al. 2006ah . 
The analysis based on a l arger sample of TeV b lazars leads 
to the same conclusion (M azin and Raue 2007T ). Recently, 
the HESS collaboration reported detection of multi-TeV 
gamma-rays from 1E S 0229+200, a BL Lac ob ject located at 
a redshift z=0.1396 (I Aharonian et al. 2007al ). It is remark- 
able that the detected hard gamma-ray spectrum of this 
source with a photon index r o b s ~ 2.5 extends up to 15 TeV. 
This, to a certain extent surprising result can be explained 
by the shape of the energy flux of EBL which between the 
NIR and MIR ba nds is expected to be proportional to A -1 
I Aharonian 20011 V Yet, the absolute EBL flux, derived from 
a rather conservative assumption that the photon index of 
the intrinsic spectrum of TeV gamma-rays does not exceed 
1.5, appears again close to the EBL lower limit, this time 



at MIR (« 2 - 3 nW /m 2 sr at 10 /mi) , derived from th e 
Spitzer galaxy counts (|Fazio et al. 20041 ; iDole et al. 20061 ). 
Thus, the gamma-ray observations of 1ES 1101-232 and 1ES 
0229+200 can be interpreted as an argument that the galax- 
ies resolved by the Hubble and Spitzer telescopes provide 
the bulk of the EBL flux from optical to mid infrared wave- 
lengths. Given the importance of such a statement, in par- 
ticular for understanding of contribution of the first stars to 
the EBL (see e.g. Kashlinsky (2005); Mapelli et al. (2006)), 
it is essential to explore alternative ways of explanation of 
very hard intrinsic gamma-ray spectra or even sharper spec- 
tral features (like pile-ups) in TeV blazars. In this context, 
recently some extreme assumptions regarding the distribu- 
tions of accelerated particles have been proposed. In par- 
ticular, Katarzynski et al (2006) argued that a gamma-ray 
spectrum as hard as To ~ 0.7 can be formed in a SSC model 
assuming a narrow parent electron distribution, e.g. power- 
law within Ei and Ei, with a low-energy cutoff Ei not much 
smaller than the high energy cutoff, Ei. In similar lines, 
Stecker et al (2007) argued that electron spectra with power 
law index ^ 1 can be accommodated within the models of 
relativistic shock acceleration . It should be noted, however, 
that in compact objects relativistic electrons usually suf- 
fer very fast synchrotron losses, therefore the assumptions 
about hard electron acceleration cannot yet guarantee hard 
gamma-ray spectra. Indeed, the radiatively cooled electron 
spectrum cannot be harder than dN/dE oc E~ 2 , indepen- 
dent of the initial (acceleration) spectrum (see e.g. Aha- 
ronian 2004). If so, the inverse Compton scattering would 
result in a gamma-ray spectrum steeper than E~ 15 . In fact, 
the Klein-Nishina effect makes the spectrum even steeper. 
In principle, one can avoid the synchrotron cooling of elec- 
trons, e.g. in a cold ultrarelativi stic wind. However suc h a 
hypothesis suggested for Mkn 501 (| Aharonian et al. 200*21 ). in 
analogy with pulsar winds, needs thorough theoretical stud- 
ies to clarify whether such cold ultrarelativistic winds can 
be formed and survived around supermassive black holes in 
the cores of AGN. 

In this paper we suggest a new scenario which allows for- 
mation of very (in practice, arbitrary) hard gamma-ray spec- 
tra in a quite natural way. The model is based on a postu- 
lation that gamma-rays before leaving the source suffer sig- 
nificant photon-photon absorption due to interactions with 
dense radiation fields inside or in the vicinity of compact 
gamma-ray production region(s). Interestingly, the presence 
of high density radiation fields of different origin in the in- 
ner parts of blazars generally is treated as a problem for the 
escape of high energy gamma-radiation from their produc- 
tion region, and, in this regard, the current models of TeV 
blazars are designed in a way to avoid the internal gamma- 
ray absorption. Below we show that, in fact, a moderate 
internal photon-photon absorption can be a clue to the very 
hard intrinsic energy spectra of TeV blazars. 



2 INTERNAL ABSORPTION OF 
GAMMA-RAYS IN BLAZARS 

When propagating through an isotropic source of low- 
frequency radiation, the gamma-ray absorption at photon- 
photon interactions is characterized by the optical depth 
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-(E) = 



a(E, e)n p h(e, r)dedr 



(1) 



where n p h(e,r) describes the spectral and spatial distribu- 
tions of target photons in the source of size R. With a 
good accuracy, the total cross-section in the monoenergetic 
isotropic photon field can be represented in the form (see 
e.g. Aharonian 2004): 



3<tt 
2s 5 " 



(2) 



The cross section depends only on the product of the pri- 
mary (E) and target photon (e) energies, s = Ee/rn^c 4 . 
Close to the threshold, s — > 1, the pair production 
cross-section behaves as <t 77 ~ (1/2)<tt(s - 1) 3/2 . The 
cross-section decreases with s also when s 3> 1: <t 77 ~ 
(2/3)cr T s _1 lns. The cross-section achieves its maximum at 



3.5: <t 7 



0.2cr T . 



For a homogeneous source with a narrow spectral distri- 
bution of photons, for order of magnitude estimates one can 
use the approximation t(E) ~ Ra 1 -y(E,e)n(e). In this case 
we should expect maximum absorption effect at gamma-ray 



energy E* 



3 c /e. Both at lower and higher energies, 



the source becomes more transparent, thus we should ex- 
pect a quite strong deformation of the primary spectrum. 
In the case t(E*) 1, the effect could be dramatic, given 
the exponential dependence of the absorption on the opti- 
cal depth. Note that while for a narrow spectral distribution 
of target photons the monoenergetic approximation gives a 
quite accurate estimate of the effect at E 3> E* , at low en- 
ergies, E ^ 1/41?*, this approximation implies a completely 
transparent source (i.e. r = 0) although, a non-negligible ab- 
sorption can take place also below E* . For example, because 
of interactions with the Wien tail, the absorption effect in 
the black-body radiation field cannot be disregarded even at 
very low energies, E <C m^c 4 /kT. 

In FigfJ] (upper panel) we present the gamma-ray at- 
tenuation factor, k = exp (— r) in a grey-body radiation 
field described by Planckian distribution with three dif- 
ferent temperatures T = 10 3 K, 10 4 K, and 10 5 K cal- 
culated for an optical depth fixed at the energy corre- 
sponding to the maximu m absorption, E* ~ m^c 4 /kT 
l|Gould and Schreder 19671 ): Since the optical depth is a 
function of the product E ■ kT, three curves are identical, 
but shifted relative to each other by a factor proportional to 
the radiation temperature. In FigfJ] we show also the atten- 
uation factors for a fixed temperature T = 10 4 K calculated 
for three different optical depths r max = 0.5, 3, and 6. It 
is seen that the gamma-ray attenuation, starting from the 
energy E ~ O.li?*, gradually increases up to E ~ E* , af- 
ter which the source becomes more and more transparent 
(k oc exp [— E~ x InE] — > 1), and, consequently, the primary 
spectrum starts to recover. As a result, in this energy inter- 
val the spectrum appears harder than the primary spectrum. 
The bottom panel of Fig[T] shows the change of the slope of 
gamma-ray spectrum, AF, which can be interpreted as a 
change of the local photon index, assuming that the initial 
gamma-ray spectrum is described by a power-law distribu- 
tion, dN/dE oc E' r °. The emerging spectrum of gamma- 
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Figure 1. Upper pannel: Attenuation factor k = exp (— r) for 
three different temperatures of target photons, T = 10 3 K, 10 4 K, 
and 10 5 K (for the optical depth r max = 6). For T = 10 4 K cal- 
culations are performed for three different optical depths: T max = 
0.5 (blue), 3 (green), 6 (red). Bottom pannel: Variation of the 
local photon index of the gamma-ray spectrum. 



rays in the energy interval (0.1 — 1)E* is steeper than the 
initial spectrum (Ar ^ 0); it recovers at E — E* (AF = 0), 
and at energies E ^ E* the spectrum becomes significantly 
harder than the initial spectrum (Ar ^ 0). For example, in 
the case of initial gamma-ray spectrum E~ 2 and r max = 6, 
the emerging spectrum of gamma-rays in the energy inter- 
val (1 - 10)£* can be very hard with T = r + Ar ~ 0, 
although the absolute flux is suppressed by an almost three 
orders of magnitude at E = E* , and an order of magnitude 
at E = 10-B*. 

Note that the requirement of a narrow spectral distri- 
bution of the target photons is a key condition for this re- 
markable effect. It should not necessarily be a Planckian or 
monoenergetic distribution, but may have any other shape, 
for example power-law with a low-energy cutoff: n(e) oc E~ a 
at e ^ Ei, and n(e) = at e ^ e\. In this case, the low- 
energy cutoff £1 plays a similar role as the temperature in 
the Planckian distribution. This is demonstrated in Fig. [2] 
for a power-law distribution of the background field with 
a = 2 and sharp cutoff at £1 = 1 eV and 10 -3 eV. Indeed, 
for the same r max = 6, the case of £1 = 1 eV is quite sim- 
ilar to the case of Planckian distribution with temperature 
T = 10 4 K shown in Fig[T] The main difference appears 
in the low-energy part, E E* . The absorption curve in 
Figf2] at low energies is smoother, because the n(e) oc £~ 2 
type distribution provides more high-energy target photons 
compared to the Wien tail of the thermal distribution, for 
interactions with low-energy gamma-rays. 

The hardening of the initial spectrum caused by inter- 
nal absorption compensates, to a large extent, the steepen- 
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Figure 2. Upper pannel: Attenuation factor calculated for a 
power-law distribution of target photons with a low-energy cutoff: 
n p h oc e~ 2 for e\ < e < oo and n p h = for e s\. Solid curves: 
£1 = 1 eV; dashed curves: e\ = 10~ 3 eV. Three different optical 
depth are shown: T max = 0.5 (blue), 3 (green), 6 (red). Bottom 
pannel: Variations of the local photon index of the gamma-ray 
spectrum. 

ing of the spectrum due to intergalactic absorption. This is 
demonstrated in Figs(3]and|4] For the EBL flux we us a "ref- 
erence" shape close to the one calculated by Primack et al. 
(2005), but with two different absolute flux normalizations 
at the wavelength A = 2.2/im: mebl(2.2 /im) = 16 nW/m 2 sr 
(Fig© and 32 nW/m 2 sr (FiggJ). The first flux is a fac- 
tor of two larger than the low-limit of EBL correspond- 
ing to the integrated light co ntributed by resolved galax- 
ies (|Madau and Pozzetti 20001 ) . while the second flux can 
be treated as an upper limit at 2.2/im (it is slightly higher 
than the fluxes c laimed from the CQBE/DIRBE and 2MAS S 
measurements (IWright et al. 200ll ; ICambresv et al. 200lf ). 
Note that for the primary (unabsorbed) differential gamma- 
ray spectrum with a photon index To = 2, the attenuation 
factor k — exp (— r) describes the SED of the absorbed ra- 
diation {E 2 dN/dE oc k(E)). All curves are obtained for 
a source located at z = 0.186. This is the redshift of the 
BL Lac object 1ES 101-232, the gamma-ray observations of 
which by the HESS collaboration have been initially used 
to constrain the EBL fl ux at optical and NIR wavelengths 
| Aharonian et al. 2006al ). 

The solid curves in Figs(3] and [4] (marked as "d") corre- 
spond to the pure intergalactic effect (i.e. without the inter- 
nal absorption). They show strong steepening of the spec- 
trum below 1 TeV and above 10 TeV with a noticeable recov- 
ery of the initial shape around a few TeV, which is explained 
by the specific shape of the EBL energy flux ( close to uebl oc 
A -1 ) between 2 and 10 /im (| Aharonian 200 ll) . It is seen that 
for the absolute flux of EBL with a normalization at 2.2/im, 



mebl(2.2 fim) — 16 nW/m 2 sr, the photon index between 
100 GeV and several TeV is changes by AF ~ 2, thus, the 
intrinsic (source) spectrum should be very hard with a slope 
T = r obs - AT ~ 1, where r ob ? = 2.88 ±0.17 is the ob served 
photons index of 1ES 101-232 (| Aharonian et al. 2006ah . Pos- 
tulating that the photon index allowed by conventional mod- 
els of gamma-ray production in blazars should not exceed 
To = 1.5, an upper limit on the EBL flux, at the level 
of UEBL (2.2 fim) w 10 nW/m 2 sr has been derived by the 
HESS collaboration (Ah aronian et al. 20 06a). It is seen from 
Fig[3] that this upper limit can be readily increased by a 
factor of 1.5, alowing a substantial internal absorption of 
gamma-rays. Indeed, the "joint operation" of internal and 
intergalactic absorptions results in the changes of the slope 
of the initial gamma-ray spectrum in the relevant energy 
band by Ar ~ to 1.5 for r max = 3, and Ar from -0.5 to 
+0.5 for T max = 6. In the latter case, the overall change of 
the shape of the initial spectrum from 100 GeV to 5 TeV is 
quite small, so hard TeV gamma-ray spectra can in princi- 
ple be detected also from distant blazars. This assumption 
can solve, to a certain extent, the problem related to the 
spectra of accelerated (parent) particles, and thus (unfortu- 
nately!) relax the constraints on the EBL that can be de- 
rived from gamma-ray observations. Formally, a detection 
of not-very-steep TeV gamma-ray spectra with r o b s S= 4 
from distant (z 0.15) blazars cannot be excluded even 
for an EBL flux close to the claimed high EBL fluxes de- 
rived from the CQBE /DIRBE and 2MASS measurements 
jCambresv et al. 200 ll ). u E BL(2.2/im) w 28 nW/m 2 sr. This 
is demonstrated in FigU] At energies above 300 GeV the 
absorption of gamma-rays in EBL only increases the pho- 
ton index by Ar ^ 4, while an additional internal absorp- 
tion with T max = 10 results in AF 2. Approximately the 
same result is expected for low EBL (e.g. at the level of 
^EBL (2.2/im) = 10 nW/m 2 sr), but for sources located at 
z ^> 0.1. In this regard, the recent claim of detection of 
VHE gamma- rays from 3C 279 (z = 0.538) by the MAGIC 
collaboration (|Teshima et al. 20071 ) can be an indication of 
significant gamma-ray absorption inside the source. Other- 
wise, even for the minimum possible EBL flux, the pure in- 
tergalactic absorption would result in an extremely steep 
VHE gamma-ray spectrum with a photon index > 5. It 
is intersting to note that the UV photons of the broad- 
line emission region of a size R ~ 1 17 cm and luminosity 
L~ 2 x 10 44 erg/s (|Pian et al. 2005h can serve as a perfect 
target for internal absorption of high energy gamma-rays in 
3C 279. 



The internal absorption of gamma-rays significantly in- 
creases the energy requirements to the source. For example 
in the case of r max = 5, the internal absorption leads to 
the reduction of the observed flux at 1 TeV by an addi- 
tional factor of 10. This, combined with the intergalactic 
absorption, implies 3 orders of magnitude attenuation of 
the primary radiation (see Fig[3}. For the quiescent state 
of 1ES 101-232, the corresponding apparent gamma-ray lu- 
minosity around 1 TeV is estimated Z/TeV — 10 44 k _1 erg/s. 
For the attenuation factor k(1 TeV) ~ 10 -3 , it becomes 
huge, 10 47 erg/s, and perhaps one or two orders of mag- 
nitude even larger in the flaring states (in analogy with 
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E, eV 

Figure 3. Internal and intergalactic absorption of gamma-rays. 
Upper pannel: Attenuation factors. The internal absorption is cal- 
culated for a Planckian distribution of target radiation with tem- 
perature T = 5 X 10 4 K. Three dashed curves correspond to the 
internal optical depths T max = 6 (a) 3 (b), 0.5 (c), (d). The 
corresponding solid curves include both the internal and inter- 
galactic absorption. The intergalactic absorption is calculated for 
a source at z = 0.186 (the redshift of the BL Lac object 1ES 101- 
232), assuming a reference shape of the EBL spectrum close to 
the one calculated by Primack et al. (2005), and normalized to 
the EBL flux at 2.2 /im: mebl(2.2 /im) = 16 nW/m 2 sr. Bottom 
pannel: Variation of the local photon index. 

Mkn 421, Mkn 501 and PKS 2155-304)0. Nevertheless, since 
there is little doubt that gamma-rays are produced in rel- 
ativistically moving jets with a Doppler factor S ^ 30 or 
even ^ 100, as it f ollows from the rec ently reported vari- 
ability of Mkn 501 (lAlbert et al. 20071 ) and PKS 2155-304 
l|Aharonian et al. 2007bl ) on minute scales (see e.g. Fabian et 
al. (2007)), the intrinsic luminosity L ln t = Z/ app <5 -4 , could be 
quite modest, namely, at the level of 10 39 — 10 41 erg/s which, 
in fact, is not far from the TeV ga mma-ray luminosity of th e 
nearby non-blazar type AGN M87 (|Aharonian et al. 2006bl ). 

The attenuation of TeV gamma-rays dramatically in- 
creases the intrinsic TeV to GeV gamma-ray flux ratio. This 
does not, however, contradict the GeV flux upper limits 

1 The detected gamma-ray flux around 200 GeV is an order of 
magnitude larger than at 1 TeV, however, because of dramatic 
reduction of the absorption effect, the contribution of low energies 
to the (absorption-corrected) apparent luminosity, is relatively 
small. 




E, eV 

Figure 4. The same as in Fig|3] but the curves are calculated 
for the internal optical depths 10 (a), 5 (b), 1 (c) and (d). For 
EBL flux is assumed twice larger than in Fig|3] «ebl(2.2 /im) = 
32 nW/m 2 sr. 



available from the EGRET observations (typically, at the 
level of 10 -10 erg/cm 2 s or higher), especially if one takes 
into account that the energy spectra of gamma-rays pro- 
duced in some principal radiation processes (e.g. through in- 
verse Compton scattering or proton synchrotron radiation) 
at sub- TeV energies are expected harder than E~ 2 . On the 
other hand, we certainly expect GeV gamma-rays from TeV 
blazars, and in this respect, the upcoming GLAST measure- 
ments with significantly improved (compared to EGRET) 
sensitivity, especially at multi-GeV energies, should provide 
the first effective probes of TeV blazars in the MeV/GeV 
domain in general, and for the "internal gamma-ray absorp- 
tion" scenario, in particular. 

Finally, in the case of hadronic origin of TeV gamma- 
rays produced at pp and/or py interactions, the flux of ac- 
companying TeV neutrinos, which freely penetrate through 
the internal and extragalactic radiation fields, can be as high 
as 10 -10 neutrinos/cm 2 s, i.e. above the detection threshold 
of the next generation km 3 scale neutrino detectors. The de- 
tection of both gamma-rays and neutrinos from TeV blazars, 
and the comparison of fluxes of these two components of ra- 
diation would provide principal information about the high 
energy processes in blazars, as well as about the attenuation 
of gamma-rays due to the (combined) internal and inter- 
galactic absorption. An additional information about the in- 
ternal photon-photon absorption alone (separated from the 
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Figure 5. Gamma-ray spectra caused by the internal photon- 
photon absorption (dotted curves) or formed during a develop- 
ment of pair cascades (solid curves) in a radiation field with tem- 
perature T = 10 4 and optical depth T max = 6. It is assumed 
that the energy density of radiation significantly exceeds the en- 
ergy density of the magnetic field, u T ^> ub- The spectra of pri- 
mary gamma-radiation (dashed lines) are assumed as "power- 
law with exponential cutoff": AN/AE oc E~ 2 exp (— E/Eo) with 
E = 10 TeV and E = oo. 

extragalactic absorption) is contained in the radiation of sec- 
ondary (pair-produced) electrons. 



3 RADIATION OF SECONDARY ELECTRONS 

The propagation of high energy gamma-rays through a low- 
energy photon field cannot be reduced to the simple effect 
of absorption. When the gamma-ray photon is absorbed, its 
energy is transfered to the electron-positron pair. The sec- 
ondary electrons interacting with the ambient magnetic and 
radiation fields produce high energy photons, either via syn- 
chrotron radiation or inverse Compton scattering. The syn- 
chotron photons are produced with much smaller energies, 
thus they do not interact with the background low-energy 
photons. In a target field with narrow spectral distribution, 
inverse Compton scattering of the photo-produced electrons 
proceeds in the Klein-Nishina limit; the upscattered photon 
receives the major fraction of the electron energy, thus is 
able to interact again with background photons. The second 
generation pairs again produce gamma-rays, thus an elec- 
tromagnetic cascade developes. 

While the energy of gamma-rays interacting with EBL 
dissipates in the intergalactic medium, and in this way con- 
tributes to the diffuse extragalactic background radiation, 
the secondary radiation caused by internal absorption may 
accompany the primary (unabsorbed) fraction of gamma- 
rays. In this regard, the development of an electromagnetic 
cascade is not a desirable process, because it masks the dis- 
tinct absorption features, and thus prevents the formation 
of very hard gamma-ray spectra. This is demonstrated in 
Fig[5] The cascade spectra are calculated assuming that the 
region of production of primary gamma-rays is located in the 
center of a spherical source filled with grey-body radiation 
with temperature T — 10 4 K and optical depth r max = 6. 
The spectrum of primary gamma-rays is given in the form: 
AN/dE oc E~ 2 e~ E/E ° for two values of the cutoff energy: 
E Q = 10 TeV and E = oo. In Fig. \5\ both the absorbed 
and cascade gamma-ray spectra are shown. It is seen that 



the development of pair cascades fully washes out the ab- 
sorption features, and instead forms a standard spectrum 
with a maximum ("bump") around the interaction thresh- 
old, ~ (m e c 2 ) 2 / kT ~ 100 GeV, followed by a steep spectrum 
above the "bump" . The cascade spectrum around 1 TeV sat- 
urates at the level of 10 percent of the primary gamma-ray 
flux. The spectrum above 1 TeV has a flat shape until the 
efficiency of the cascade drops (because of the progressively 
decreasing optical depth) with a gradual transition to the 
"absorption" regime. 

There are two ways of significant reduction of the contri- 
bution from the cascade component to the observed gamma- 
radiation. 

(i) Absorption of gamma-rays outside the production re- 
gion. In this case the "foreground" cascade radiation cannot 
screen the unabsorbed fraction of gamma-rays. Indeed, al- 
though the energy in the cascade radiation exceeds, by a fac- 
tor of 100, the energy of the unabsorbed fraction of primary 
radiation (see Fig. [5}, for the observer the primary radiation 
emitted by a relativistically moving source ('blob") will be 
much brighter (oc t5|) and shifted to higher energies (oc Sj) 
compared to the isotropic cascade emission of the surround- 
ing environment. On the other hand, the part of the cascade 
developed inside the blob is, of course, also Doppler boosted, 
therefore the condition of suppression of the cascade com- 
ponent can be satisfied when the optical depth within the 
blob t' <C 1, i.e. the gamma-ray production region is much 
smaller than the source of the optical/UV radiation, I <C R- 
Since th gamma-ray production regions in blazars are be- 
lieved to be very compact, I ~ 10 14 — 10 1B cm, we may 
conclude that the source of the optical radiation should be 
larger than R ~ 10 15 — 10 17 cm. There are many poten- 
tial sources of optical radiation in blazars (see e.g. Urry and 
Padovani 1995). In this paper we do not intend to spec- 
ify the origin of the low-energy radiation fields, but simply 
notice that an even modest optical source located in the 
core of a blazar can provide an effective target for gamma- 
ray absorption. The luminosity of this source is estimated 
as Lo = l2irn p i l kTR 2 c, where n p h is the average number 
density of radiation. It can be estimated from the condi- 
tion a 77 i?n p h = T max . For a characteristic optical depth 
Tmax ~ 5 and temperature T — 5 x 10 4 K, we obtain, 
Lo ~ 2 x 10 43 (i?/10 17 cm) erg/s. 

(ii) Secondary electrons cooled through synchrotron radia- 
tion. This condition can be satisfied if the energy density 
of the magnetic field exceeds the energy density of radia- 
tion, B 2 /8ty > 3kTn ph or B > 0.4(i?/10 17 cm)" 1/2 G. In 
Fig. [6] the broad-band SED of the radiation initiated by 
absorption of primary gamma-rays with a power-law spec- 
trum AN/AE oc E~ 2 is shown, assuming that the absorption 
of gamma-rays takes place inside the gamma-ray production 
region. Calculations are performed for three different optical 
depths T ma x=0.5, 3, and 6, a fixed temperature of radiation 
T = 5 ■ 10 4 K, and two extreme values of the magnetic field, 
B — 100 G and B — 0.1 G. In calculations of the electron 
spectra we assume that the energy losses of electrons are 
dominated by synchrotron cooling. Note that the spectrum 
of synchrotron radiation of secondary electrons has a char- 
acteristic bell-type form, i.e. quite similar to the synchrotron 
spectra formed in the synchrotron-self Compton (SSC) mod- 
els. However, in the SSC models designed for TeV blazars 
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Figure 6. Synchrotron radiation of secondary electrons. The pri- 
mary gamma-ray spectrum is assumed dN/dE oc E~ 2 (dotted 
line). The target photon field is Planckian with T = 5 ■ 10 4 K. 
The absorbed gamma-ray spectra and the corresponding spec- 
tra of secondary electrons are calculated for three optical depths 
T max = 6 ( a )i 3 (b), and 0.5 (c). The synchrotron radiation is cal- 
culated assuming that the absorption of gamma-rays takes place 
inside the source (gamma-ray production region) for two values 
of the magnetic field: B = 100 G (solid curves) and B = 0.1 G 
(dashed curves). 

the synchrotron peak is a result of the maximum energy of 
accelerated electrons, while the hard low-frequency part of 
the spectrum is determined by radiation of uncooled low- 
energy electrons. In the "internal gamma-ray absorption" 
scenario we see similar features, but for different reasons. 
The electrons in a narrow-band radiation field are produced 
with a spectrum similar to the spectrum of parent gamma- 
rays (oc E~ T °), but with cutoffs both at low- and high en- 
ergies. These cutoffs are explained by the threshold of the 
photon-photon pair production and by the reduction of its 
cross-section at highest energies, respectively. Due to the 
synchrotron cooling, at low energies the electron spectrum 
obtains a standard oc E~ form, which gradually transforms 
to a oc E~ ( - r ° +1 ^ type spectrum at intermediate energies and 
a cutoff at highest energies. The corresponding spectral fea- 
tures are reflected in the synchrotron spectrum - a power- 
law with a photon index 1.5 at low energies, with a smooth 
transition to a power- law with a photon index To /2 4-1 at in- 
termediate energies, and a smooth gradual cutoff at highest 
energies. A significant difference between the SSC and the 
"internal gamma-ray absorption" models appears also in the 
ratio of fluxes corresponding to the low (IR to X-ray) and 
high (gamma) energy peaks. While in the SSC scenario this 
ratio is determined by the ratio of energy densities of the 
magnetic and target photon fields, in the "internal gamma- 
ray absorption" scenario this ratio is basically determined 
by the efficiency of gamma-ray absorption inside the source. 
The magnetic field determines only the position of the syn- 
chrotron peak, but not the flux level. The latter is deter- 
mined by the factor proportional to (1 — exp [— T max ]). For 
Tmax ^ 1, the dependence on the optical depth almost dis- 
appears. All these features can be seen in Fig. [6] 

The position of the synchrotron peak depends, in a quite 
interesting way, also on the Lorenz and Doppler factors of 
the source. While in the standard models of blazars the 
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Figure 7. Impact of the bulk motion on synchrotron radiation 
of secondary electrons. The magnetized blob moves with a rela- 
tivistic velocity. The synchrotron maximum moves by a factor of 
5j'/7? and the distributions become somewhat wider. The pri- 
mary gamma-ray spectrum was assumed power-law: oc E7 2 (is 
shown by the dotted line). The target photon field is Planckian 
with T = 5 ■ 10 4 K. The absorbed spectra are shown by black 
solid curves for r max = 6(a), T max = 3(b) and r max = 0.5(c). We 
calculate synchrotron radiation from a magnetized region with 
the average optical depth (over different directions) = 6(a, 

red lines), 3(b, green lines) and 0.5(c, blue lines). The magnetic 
field is assumed B = 100 G. Calculations are performed for two 
different bulk Lorenz factors: jj = l(no motion) and jj = 30 (for 
5j = 20). 

Doppler effect shifts the overall SED towards higher ener- 
gies by a factor of Sj, in the "gamma-ray absorption" sce- 
nario the synchrotron peak is shifted towards lower ener- 
gies (see Fig[7| The reason is quite simple. In the frame of 
the relativistically moving source illuminated by an external 
radiation with a characteristic energy eo, the electrons are 
produced with energies E min ^ (m e c 2 ) 2 /(e 7j), thus the po- 
sition of the synchrotron peck in the frame of the observer 
is proportional to SjE^in oc 5/y 2 . Thus in a source moving 
towards the observer at a small angle, the position of the 
synchrotron peak in the frame of the observer is inversely 
proportional to the Doppler factor 5j, just opposite to the 
spectrum of gamma-rays which is shifted towards higher en- 
ergies by the same Doppler factor. 

The position of the synchrotron peak depends strongly 
also on the average energy (or temperature) of the target 
radiation field. Indeed, with an increase of the temperature 
of background radiation T, the threshold of photon-photon 
interactions, and consequently the minimum energy of pro- 
duced secondary electrons decreases as i? m in oc 1/T, and 
hence the synchrotron peak moves towards lower energies 
as hvm oc 1/-E m i n oc 1/T 2 . Generally, the position of the 
synchrotron peak of pair produced electrons depends on the 
temperature of the target radiation field, the magnetic field, 
and the Doppler and Lorentz factors of the jet, as 

hu m oc BT- 2 {5 3 /^) . (3) 

It is easy to derive a simple analytical expression which de- 
scribed the high energy part of the synchrotron spectrum, 
hv 2> hvm, produced by electrons for which the source be- 
comes optically thin, t(E) ^ 1. In this case the produc- 
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tion spectrum of electrons Q(E) oc 1/ EAN 1 /AE (here we 
ignore the weak logarithmic term in the photon-photon in- 
teraction cross-section). Then, for the power-law gamma-ray 
spectrum, AN-y/AE oc E~ F °, the cooled electron spectrum 
is also power-law, AN e /AE e oc E~ r "~ 2 , and correspond- 
ingly the SED of synchrotron radiation, vF v oc ^"To+i-)/ 2 ^ 
For example, for To = 2, the SED of synchrotron radia- 
tion is rather steep, vF v oc v~ . In fact, because of the 
cutoff in the gamma-ray spectrum, the high energy tail 
of synchrotron radiation is expected even steeper. This is 
demonstrated in Fig. [8] where the broad-band SEDs of ra- 
diation initiated by gamma-rays in a source at z — 0.186 
are shown. It is assumed that gamma-rays in the frame 
of the jet moving with a Lorentz factor jj — 10 have a 
power-law distribution with an exponential cutoff at 1 TeV, 
diV 7 /d£; oc £'" 3/2 exp(-£/lTeV). It is assumed also that 
Sj — 7j (i.e. the viewing angle is 9 w 6°). The calcula- 
tions are performed for two temperatures of the radiation 
field through which the jet propagates - T = 5 x 10 4 K and 
T = 5 x 10 K, assuming that in both cases the optical 
depth inside the moving gamma-ray production region (the 
blob) with a homogeneous magnetic field is r max = 3. Finally 
for the intergalactic absorption a "template" EBL spectrum 
is assumed with a normalization at 2.2 /im at the level of 
16 nW/m 2 ster. The impact of the temperature on both the 
spectrum of arriving gamma-rays and the synchrotron radi- 
ation of secondary electrons is clearly seen. Note that while 
below 100 GeV the deformation of the primary gamma-ray 
spectrum is caused mainly by internal absorption, the sharp 
cutoff at energies above 10 TeV is due to the severe inter- 
galactic absorption. In the intermediate energy range be- 
tween 100 GeV and 10 TeV the internal and intergalactic 
photon-photon interactions "operate" together resulting in 
quite specific broad-band SEDs. The discussion of impli- 
cations of these results for specific astrophysical objects is 
beyond the scope of this paper. We note only that our pre- 
liminary studies show that the suggested model in general 
can satisfactorily explain the observed broad-band SEDs of 
TeV blazars. 



4 DISCUSSION 

The energy spectra of VHE gamma-rays from blazars, af- 
ter correction for intergalactic absorption, generally appear 
very hard, even for the minimum flux level of EBL deter- 
mined by the integrated light of resolved galaxies at optical 
(Hubble) and infrared (Spitzer) wavelengths. A slight devia- 
tion from the robust lower limits of EBL leads to unusually 
hard intrinsic gamma-ray spectra which cannot be easily 
explained within the standard particle acceleration and ra- 
diation models. In this paper we suggest a scenario which 
can lead to the formation of instrinsic gamma-ray spectra 
of arbitrary hardness without introducing modifications in 
the particle acceleration models. The main idea is that the 
gamma-rays before they leave the source suffer significant 
internal energy-dependent absorption due to interactions 
with the ambient low-frequency photons. The existence of 
dense radiation fields of different origin in blazars (see e.g. 
Urry and Padovani (1995)) combined with the large photon- 
photon pair production cross-section, makes this scenario 
quite natural and effective, in particular in the compact 



cores of blazars. For the formation of hard VHE gamma- 
ray spectra, the target radiation field must have a rather 
narrow spectral distribution or a sharp low-energy cutoff, 
with a typical energy of photons of about 1 to 10 eV. For- 
mally, for very large optical depths, this process can provide 
an arbitrary hardness of gamma-ray spectra, though at the 
expense of a significant increase of the required nonther- 
mal energy budget. However, as long as the current blazar 
models require relativistically-moving gamma-ray produc- 
tion regions with large Doppler factors, 8j ^ 30, and per haps 
even more (lAharonian et al. 2007bl ; iFabian et al. 20071 ') , the 
available energy budget seems to be not a critical issue. 

The unavoidable feature of the proposed model is the 
radiation of secondary electrons via synchrotron or inverse 
Compton scattering. If the optical depth inside the gamma- 
ray production region is small, r -C 1, e.g. the gamma-ray 
source is much smaller than the external source of optical 
photons, the secondary electrons are produced and radiate 
mainly outside the gamma-ray production region. Even in 
the case of heavy absorption of gamma-rays, the secondary 
radiation of secondary electrons can hardly be detected. 
Indeed, since the intrinsic gamma-ray luminosity is rela- 
tively modest, and the absorbed energy is re-radiated as 
an isotropic sourc^E the lost of the beaming factor dramat- 
ically reduces the signal compared to the primary (Doppler 
boosted) radiation. 

The picure is dramatically changed when the gamma- 
ray source moves through a very dense photon field, such 
that the optical depth inside the source becomes larger than 
1. In this case the main fraction of the absorbed energy is re- 
leased in the form of secondary electrons inside the gamma- 
ray production region, and thus the radiation of the sec- 
ondary electrons profits, as the primary gamma-radiation 
does, from the Doppler boosting. The secondary electrons 
are cooled through synchrotron and/or inverse Compton 
channels. The latter in fact proceeds via development of 
pair cascades as long as the typical energies of electrons or 
gamma-rays and the energy of target photons eE Cil ;§> m 2 c 4 . 
The cascade, however, diminishes the energy-dependent ab- 
sorption features, thus the model becomes effective when 
the electrons are cooled predominantly via synchrotron ra- 
diation, i.e. B 2 /87r u T . The energy density of the radiation 
u T = £7i p h with an average energy of target photons of about 
e ~ 1 eV is estimated from the condition r max ^ 1, thus for 
the effective suppression of the cascade 

B > (407re/(7T«) 1/2 « Q.5(l/10 15 cmy 1/2 T^ x G (4) 

For an optical depth r ma x ~ 1 and the size of the 
gamma-ray source I ~ 10 16 cm, the magnetic field exceeding 
0.1 G should be sufficient to prevent the cascade. For smaller 
gamma-ray production regions, e.g. I ~ 10 14 cm, the mag- 
netic field should be larger than 1 G. For such magnetic 
fields the synchrotron radiation of secondary electrons ap- 
pears in the optical to hard X-ray energy bands. Depending 
on the optical depth, the synchrotron peak can be higher 
than the gamma-ray peak. Interestingly, unlike the classi- 
cal "synchrotron/Inverse-Compton" models, where the ra- 

2 Unless the electrons are produced in an environment with a 
very low magnetic field, and thus are cooled via inverse Compton 
scattering before any noticeable deflection. 
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Figure 8. Impact of the radiation temperature on the spectral energy distribution of emission initiated by primary gamma-rays in a jet 
moving through a homogeneous radiation field. The source is located at z = 0.186. The dashed curve is the assumed primary spectrun of 
gamma-rays. The solid curves repreesent the synchrotron raduiation of the secondary electrons as well as the gamma-ray spectra after the 
internal and extragalactic absorption. Calculations are performed for two radiation temperatures, T = 5 X 10 4 K (a) and T = 5 X 10 5 K 
(b). For both cases 8j = ■yj = 10, B = 1 G, T m ax = 3. An EBL is assumed with a normalization of the flux «EBL(2.2^tm) = 16 nW/m 2 str. 



tio of the synchrotron to IC peak is determined by the ra- 
tio ub/u t , in the "internal gamma-ray absorption" scenario 
the synchrotron peak does not strongly depend on the mag- 
netic field. Whether this scenario can be applied to the the 
broad-band SEDs of gamma-ray blazars, is an interesting 
issue which requires special dedicated studies. 

Finally, we want to discuss briefly the radiation mech- 
anisms of primary gamma-radiation. Generally, the model 
does not give a preference to the leptonic or hadronic origin 
of radiation, unless the magnetic field exceeds the estimate 
given by Eq.Q. In this case the synchrotron-to-IC flux ra- 
tio produced by directly accelerated electrons would be too 
high, especially after the internal absorption of gamma-rays, 
contrary to the detected SEDs of most of the TeV blazars. 

Large magnetic fields in the gamma-ray production re- 
gion, typically B ^ 1 G, would favor gamma-ray production 
by relativistic protons, with all advantages and disadvan- 
tages common for hadronic models. The basic problem of 
hadronic models is linked to the low interaction rates which 
do not allow the most natural explanation of the observed 
fast gamma-ray variability of blazars in terms of radiative 
cooling. For example, in the case of interactions of protons 
with the ambient plasma with number density n, the charac- 
teristic time of pp interactions with production of 7r°-mesons 
is t pp « 10 15 n _1 s. Thus, in order to explain the variabil- 
ity of gamma-rays as short as several minutes like the TeV 
flares observed from PKS 21555-301 and Mkn 501, the den- 
sity of plasma should be as large as 5 x lO 12 ^" 1 cm -3 which 
implies a very heavy source and correspondingly huge ki- 
netic energy E^ u = (4/3)irl 3 nm p c 2r Yj « 10 55 erg (here we 
assume that S~yj). One may invoke alternative explanations 
of the variability of blazars, e.g. due to the adiabatic losses 
or escape of particles from the source, but this assumption 
leads to dramatic reduction of radiation efficiency, and to an 
increase the energy requirements to the accelerated protons. 

A similar problem face the photomeson processes at 
interactions of protons with the ambient radiation fields. 
Actually in the "internal gamma-ray absorption" scenario 
this mechanisms seems a quite natural choice because the 
same background photons which absorb gamma-rays can 



play a role of the target for photomeson interactions. 
However, because of the small cross-section, the efficiency 
of this process again appears quite low. The interaction 
time of protons with energy, E ^ 200 MeV/(e7j) ~ 
2 x 10 16 (e/l eV)" 1 (7 j /10)" 1 eV (in the frame of the 
moving source with a Lorentz factor jj) is estimated 

*P7 ~ l/(.M>7™ph c ) ~ (cr 77 /(7p 7 )/ _1 i?/cr m ax (< °" P 7 >~ 

10 -28 cm 2 is the average cross-section and / ~ 0.2 is the 
multiplicity of the process). Thus, we can see that during the 
passage of the source of optical photons of size R, the pro- 
tons transfer only er P7 /<7 77 ~ 10~ 3 fraction of their energy to 
gamma-rays. If such a low efficiency can be compensated by 
very large Doppler boosting (e.g. assuming 5j ~ 100), this 
channel can provide very large fluxes of neutrinos, which 
unlike gamma-rays do not suffer internal and extragalactic 
absorption. In the case of attenuation of VHE gamma-ray 
fluxes by 2 to 3 orders of magnitude, the expected fluxes 
of neutrinos from TeV blazars can be as large as the de- 
tection threshold of the km 3 volume high energy neutrino 
telescopes, F u (Js 1 TeV) « 10 -11 neutrinos/cm 2 s. 

It is interesting to note that, because of the threshold 
of photomeson production, the interactions of protons of ar- 
bitrary distribution with a narrow band radiation with a 
characteristic energy e, result in a differential gamma-ray 
spectrum which below the energy « 10 16 (e/l TeV) -1 eV is 
extremely hard, &N/&E = const, thus this process itself can 
provide very hard gamma-ray spectra independent of the 
spectrum of parent protons. 

Despite certain attractive features, this mechanism 
faces the same problem as pp interactions - a low radia- 
tion efficiency. Therefore it can work only under conditions 
of extremely large Doppler boosting of radiation. The effi- 
ciency of VHE gamma-ray production can be much higher 
in the case of synchrotron radiation of protons, provided 
that the acceleration of protons proceeds at a rate close to 
the fundamental limit, and the magnetic field in the proton 
accelerator well exceeds 10 G. In particular, in the mag- 
netic field of order 100 G, protons can be accelerated to 
energies 10 20 TeV and thus can produce VHE synchrotron 
gamma-rays on timescales of 10 4 s. Although due to the 
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self-regulated synchrotron cutoff (|Aharonian 20001 ) the spec- 
trum of gamma-rays is limited by sub-TeV energies, an ob- 
server detects Doppler boosted gamma-radiation extending 
to multi-TeV energies. The characteristic feature of this 
mechanism is the very large electromagnetic energy con- 
tained in the blob, l/6l 3 B 2 « 2 x 10 48 erg, hard X-ray emis- 
sion of the secondary (pair-produced) electrons, and negli- 
gible fluxes of neutrinos. 
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